Making use of geometrical invariants in black hole collisions. 
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We consider curvature invariants in the context of black hole collision simulations. In particular, we 
propose a simple and elegant combination of the Weyl invariants / and J, the speciality index S. In 
the context of black hole perturbations S provides a measure of the size of the distortions from an 
ideal Kerr black hole spacetime. Explicit calculations in well-known examples of axisymmetric black 
hole collisions demonstrate that this quantity may serve as a useful tool for predicting in which cases 
perturbative dynamics provide an accurate estimate of the radiation waveform and energy. This 
makes S particularly suited to studying the transition from nonlinear to linear dynamics and for 
invariant interpretation of numerical results. 
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I. INTRODUCTION 

Einstein's theory of gravity demands the equivalence 
of all coordinate representations of gravitational dynam- 
ics. This coordinate gauge invariance makes general rel- 
ativity very simple and beautiful because there are no 
special families of observers to be considered. On the 
other hand, though, it can also make it difficult to distin- 
guish whether differences observed in two spacetime rep- 
resentations are true physical (geometrical) differences 
or gauge differences. A natural way to limit confusion 
between gauge and physical differences is to work, when- 
ever possible, with geometrically defined scalars, which 
are invariant under (passive) coordinate transformations. 
Geometric curvature invariants have a long productive 
history in the classification and distinction of exact an- 
alytic solutions of Einstein's equations, particularly for 
the algebraic Petrov classification and characterization of 
curvature singularities 0. To a lesser extent, curvature 
invariants have also been applied in the field of numerical 
relativity, primarily for code testing when evolving ex- 
act solutions numerically. Here gauge invariant methods 
have the distinct advantage that they can be applied to 
evolutions using numerically generated coordinates which 
are not understood analytically. What seems to have re- 
ceived less attention is the application of curvature in- 
variants in a region of spacetime which is perturbatively 
close to an exact solution of Einstein's equations. 

Black hole perturbation theory has recently generated 
much interest as a model for the late stages of a black 
hole collision spacetime B. When two black holes are 
close to each other one can simply treat the problem as 
a single distorted black hole that 'rings down' into its fi- 
nal equilibrium state. Perturbative calculations applied 
in this final regime have become an important tool in 
the verification and interpretation of numerically gener- 
ated results IjH- More ambitiously, the perturbative 
approach can be used in conjunction with full scale 3D 
numerical relativity simulations to directly "take over" 
and continue a previously computed numerical black hole 



spacetime |6[. It is in the context of such an approach 
to black hole collisions that the authors expect to make 
special use of curvature invariants. 

In setting up such a perturbative approach one encoun- 
ters a family of Cauchy data sets which encode distorted 
black holes such as a family of slices from the late-time 
region of a numerically calculated black hole coalescence. 
We would expect, on physical grounds, that for some 
of the initial data sets, perturbation theory can provide 
an accurate model of the corresponding future- evolved 
spacetime. We would like to have a working criteria for 
when we can expect perturbation theory to be effective 
based only on numerical data. 

Motivated by this purpose we introduce an invariant 
quantity S, which geometrically measures local deviation 
from algebraic speciality, and which we expect to be a 
very useful tool for numerical and perturbative work in- 
volving near-stationary regions of black hole spacetimes. 
Other indicators of the potential success of perturbation 
theory, like for instance the size of the distortions of the 
apparent horizon, have been previously applied to nu- 
merical results 0], but these require establishing a coor- 
dinate system by an intuitively reasonable method and, 
in most of the cases, are computationally expensive and 
only applicable when perturbations are reasonably small. 
Our invariant index is simple, elegant and has none of 
these shortcomings. In fact, it is not at all limited to 
perturbation studies, but can be applied directly to full 
3D numerical evolutions to directly explore the transition 
from nonlinear to the linear regime |6|, or for invariant 
interpretation of numerical spacetimes. In all well-known 
examples of axisymmetric black hole collisions we stud- 
ied, <S has already proven to be very useful. 

II. THE SPECIALITY INDEX 

Curvature invariants are part of the standard analysis 
of exact solutions of Einstein's equations. From the Weyl 
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tensor, C a bcd, which carries information about the grav- 
itational fields in the spacetime, one can algebraically 
derive two complex curvature invariants usually called I 
and J. These are essentially the square and cube the self- 
dual part, C abcd = C abcd + («/ 2 ) £ q. 
tensor: 
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These scalars are useful in the algebraic classification 
of exact solutions. The different algebraic Petrov types 
are distinguished by the degeneracies among the (up to) 
four principal null directions (PNDs) associated point- 
wise with the Weyl tensor. Type I is the algebraically 
general case with four distinct principal null directions. 
The other types II, III, D, and N have at least two coinci- 
dent PNDs and are referred to as algebraically special. A 
notable characteristic common to all stationary isolated 
black hole solutions of general relativity is that they are 
all algebraically special, of Type D, with two pairs of co- 
incident principal null directions (PNDs) at each point. 
We contend however that for interesting cases involving 
nontrivial dynamics, perturbed black hole spacetimes are 
generically not algebraically special, of type I. 

Significantly, the invariants I and J satisfy the rela- 
tion, I 3 = 27 J 2 if and only if the Weyl tensor is alge- 
braically special Q . Since the Weyl tensor in a perturbed 
black hole spacetime is not expected to be algebraically 
special, we expect J 3 ^ 27 J 2 . Our proposal, then, is 
to use this relation to construct an invariant index for 
algebraic speciality as a local measure of the size of the 
distortions from some background black hole. This vi- 
olation can be in general quantified by considering the 
following speciality index 
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For the unperturbed algebraically special background 
Kerr spacetime 5=1. In the perturbed spacetime we 
generically expect 5=1 + AS, and the size of the devi- 
ation AS ^ can be used as a guide to predicting the 
effectiveness of black hole perturbation theory. 

The theory of perturbations on a background Kerr 
spacetime was worked out first by Teukolsky Jq] and has 
been extensively studied by many authors f9Hl2|l . In this 
context it is natural to use the Newman-Penrose decom- 
position of the Weyl tensor into five complex quantities, 
xpo, ipi, ip2, 03, an d 04, defined with respect to some 
choice of a null tetrad basis. In terms of the Weyl com- 
ponents, for an arbitrary tetrad choice: 
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For any Type D spacetime such as Kerr, a tetrad can be 
conveniently chosen such that only 02 is non- vanishing. 
By expressing S with respect to any perturbation of such 
a tetrad we find that 
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were e is a perturbation parameter, and the superscript 
(0) and (1) stand respectively for background and first 
order pieces of the perturbed Weyl scalars. Thus, the 
lowest order term in the deviation is second order in the 
perturbation parameter e. In the perturbative context 
this means that, when the speciality index is significantly 
different from unity, one can see that a potentially sec- 
ond order quantity has become significant, and the first 
order theory should no longer be trusted. For the case of 
Schwarzschild Eq. (|j) can be reexpressed in terms of the 
gauge invariant Moncricf functions. p0(]. 



III. EXAMPLE APPLICATIONS 

Since S = 1 in the background, a reasonable crite- 
rion for when we can expect perturbation theory to be 
applicable is that S should differ from one by no more 
than a "factor of two" . In regions where this condition is 
violated we can expect significant violations of the per- 
turbative dynamics. As test cases, to verify our interpre- 
tation of S we have considered some well-studied cases 
of axisymmetric black hole collisions. We find it useful 
to consider the location of such violating regions with 
respect to the background black hole horizon and the 
perturbative "potential barrier." As is well-known, per- 
turbative black hole dynamics are governed by a wave 
equation with a potential concentrated in the vicinity of 
r = 2M in the isotropic coordinates used for our exam- 
ples. The potential has the effect of largely preventing 
waves from crossing this region. 

The examples we consider here all correspond to even- 
parity modes implying that S is real, so we leave out 
reference to its imaginary component in the following 
discussion. We first consider the case of two initially 
resting equal-mass black holes. This initial configuration 
is represented by the equal-mass time-symmetric Misner 
datasets parameterized by /io as a measure of the ini- 
tial separation. In this case, Schwarzschild black hole 
perturbation theory has been shown to provide a very 
good estimate of the total radiated energy for cases with 
/io < 1.8 even though the black holes share a common ap- 
parent horizon only when fi Q < 1.36. Comparisons with 
numerical calculations [Q and second order calculations 
|l3f have demonstrated that the linear perturbation ap- 
proach overestimates the radiation energy by only a fac- 
tor two up to /io = 1.8 but beyond that the differences 
grow quickly. Second order perturbations have been ap- 
plied in this case as a useful tool for assessing the do- 
main of validity of perturbation theory. We can obtain 
similar conclusions by applying our speciality index test. 
Fig. |l| shows the initial values of S along the equator for 
Misner data, at several initial separations. In isotropic 
coordinates, as used in previous studies, the horizon of 
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FIG. 1. The equatorial values of S shown in this plot al- 
low the interpretation that perturbation theory can work for 
cases of fio < 1.8 but not higher. The "potential barrier" 
near ri ao = 2M prevents the violation of perturbative dy- 
namics inside near the horizon at r — 0.5M, from affecting 
the radiation. For larger separations the interpretation of a 
potential barrier at n so = 2M is itself invalidated. 

the background black hole is located at r — 0.5M and 
location of the perturbative "potential barrier" is near 
r = 2.0M. First consider the case, /io = 1.2. While 
there is small region near the horizon where the criterion 
fails, this region is well inside the potential barrier which 
should prevent any local errors in the perturbative dy- 
namics from having a significant effect on the outgoing 
radiation. In the marginal /i = 1.8 case, the error in S 
begins to be significant near the potential barrier, and 
the radiation should be somewhat affected. For larger 
values of /j,q the violation is significant in the vicinity of 
the potential barrier itself, invalidating the perturbation 
dynamics. The sudden drop in the value of S in these 
cases makes our interpretations insensitive to the choice 
of a "factor of two" cutoff. 

Another well-studied case is the collision of black holes 
with non-vanishing initial linear momentum P. Here 
again, first order perturbation theory has been very suc- 
cessful, unexpectedly providing a good estimate of the 
radiation energies even for large values of P [Q. We 
consider configurations corresponding to a fixed /i = 1.5, 
for various values of P. The corresponding initial values 
of S are shown in Fig. ^|. The presence of momentum in 
the initial slice introduces qualitatively different features 
to the initial values of S exhibiting now a region of iS > 1 
which falls off more slowly at large r. There is no ques- 
tion in this case of the location of this region, but rather 
the magnitude of the violation. In this case the violation 
seems to grow quadratically with P and reaches a factor 
of two just after P = M. Suggesting perturbation theory 



FIG. 2. The initial inward linear momentum of the black 
holes results in qualitatively different features in the equa- 
torial S values. These curves correspond to separations of 
fj,o = 1.5 and suggest that perturbation theory should be good 
below P = M. 

should be successful up to this P ~ M. The prediction of 
our prescription is consistent with the accuracy to 10% 
for radiation energies shown in Ref. |1J] for P < M but 
does not explain the mere doubling of this discrepancy 
out to P = 3M. 

In the above analysis we have considered the location 
of the 5^1 regions in relation to the local proper- 
ties of the "background" black hole. The specification 
of the background metric in this standard treatment is 
not itself gauge independent though. We can evaluate 
the applicability of a perturbative treatment in a gauge 
invariant way by utilizing a gauge-invariant specification 
of the background. A simple way to do this is to spec- 
ify the background Schwarzschild radial coordinate by 
r Schw = 3M 2 //. The location of the horizon and po- 
tential barrier are then found with respect to this coor- 
dinate. A two dimensional representation of the results 
for Misner initial data is given in Fig. ||. These plots 
show three curves, representing the locations of the back- 
ground horizon, the potential barrier, and the S = 0.5 
surfaces in a quadrant of the xz-plane. The qualitative 
features observed in the preceeding interpretation of S 
for the Misner problem, are reproduced precisely in this 
more complete, fully invariant perspective. 

IV. DISCUSSION 

We have identified a gauge invariant quantity which 
provides a particularly interesting local reduction of the 
geometric data in a (numerical) black hole spacetime. 
We foresee three areas of application where S should be 
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FIG. 3. Here we reexamine the Misner results, for the 
fio ~ 1-2, 1.4, 1.8,2.2 cases respectively, showing the full 2D 
picture in a quadrant of the xz plane. The results of the 
analysis this time provided in terms of an invariantly defined 
background and rschw = 2M horizon and potential barrier 
near rschw = 3M. The solid curve indicates the 5 = 1/2 
surface. 



a useful quantity: perturbation studies, which we have 
discussed in detail, numerical spacetime interpretation, 
and numerical code testing. In the context of pertur- 
bation theory we have demonstrated that S provides an 
invariant criterion for predicting when perturbation the- 
ory might provide a reliable approximation for part of 
a black hole spacetime. This method has the advantage 
over other estimates such as apparent horizon formation 
that it is genuinely gauge invariant. In practice it can 
serve as simple alternative to second order perturbation 
theory, but our prediction does not provide such a direct 
validation of a perturbative calculation. Knowing that 
5^1 does not, for example, identify the appropriate 
background spacetime, a vital step in any application, 
but it is very useful predictor of when perturbation the- 
ory may be a useful alternative to numerical simulation 
in a generic black hole spacetime. The quantity S itself is 
not restricted to perturbation studies but should also be 
useful in the general interpretation of numerical space- 
times. As an example, looking at S in numerical sim- 
ulations, Misner data evolve after a short time to qual- 
itatively resemble the black hole data with inward mo- 
mentum discussed above. After longer evolutions in this 
family S tends to approach unity in the horizon/potential 
region with evidence of radiation moving away in both di- 
rections. In this context we also note that the presence 
of an isolated horizon [ ^5[ , a recent construct of grow- 
ing theoretical interest in black hole spacetimes, implies 
locally that 5 = 1. Lastly, because S is an invariant 
with often predictable behavior, it can be very useful in 



numerical code testing. In Kerr spacetimes for example 
5 = 1 exactly in any coordinates. Also in typical cases of 
Bowen-York binary black hole data that we have looked 
at S falls off quickly toward unity away from the black 
holes. In light of its simplicity and straightforward sig- 
nificance, we expect S to become a standard, very useful 
tool for analyzing numerically generated spacetimes and 
interpreting their physical content. 
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